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Abstract 

Strong experimental evidence has indicated that tumor growth belongs to the molecular beam 
epitaxy universality class. This type of growth is characterized by the constraint of cell proliferation 
to the tumor border, and surface diffusion of cells at the growing edge. Tumor growth is thus 
conceived as a competition for space between the tumor and the host, and cell diffusion at the tumor 
border is an optimal strategy adopted for minimizing the pressure and helping tumor development. 
Two stochastic partial differential equations are introduced in this work in order to correctly model 
the physical properties of tumoral growth in (1+1) and (2+1) dimensions. The advantages of these 
models is that they reproduce the correct geometry of the tumor and are defined in terms of polar 
variables. Analysis of these models allow us to quantitatively estimate the response of the tumor 
to an unfavorable perturbation during the growth. 

PACS numbers: 87.10.+e, 68.35.Fx 
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In Nature, one can find a huge number of systems tliat develop a rougli interface in the 
process of growing. Many of them have been adequately understood by means of the use of 
some tools from fractal geometry, as scaling analysis, as well as via modelling with stochastic 
p„«.I d.faen.ia, e,ua«o„. (SPDE, and d..ete .ode. !i WMle t.e. „ do „ot 
constitute an exclusive theoretical framework for surfaces in the physical world, they can 
be applied to win a deeper understanding of many processes in biology P]- Due to the 
many possible important applications in medicine, tumor growth constitutes one of the 
most interesting subjects of study to which scaling analysis can be applied. Actually, a 
very important research on tumor growth has been recently carried out. It has been found 
a strong empirical evidence that a broad class of tumors belong to the same universality 
class: the molecular beam epitaxy (MBE) universality class 3, ^- MBE dynamics are 
characterized by a number of features which include a linear growth rate, the constraint of 
growth activity to the outer border of the tumor, and surface diffusion at the growing edge; 
all of them have been observed experimentally. Surface diffusion has been identified as an 
optimal mechanism for favoring tumor growth. The host tissue exerts pressure on tumors 
which opposes their growth, but surface diffusion drives the cells to the concavities of the 
interface keeping this pressure to a minimum. These findings suggested studying the effect 
of the immune response on the tumor, and it has been established that an enhancement of 
the immune response increases the pressure on the tumor surface, and therefore limits its 
development j^. A very important consequence of this fact is its possible application to 
improve cancer therapy, something that has been already exploited with positive results 0| . 
All these achivements underline the fundamental importance of understanding the physics 
of tumor growth. 

Before introducing any model, it is of fundamental importance to point out that the 
application of a physical model to a complex biological process implies a high simplification 
of many of its features. The movement of the cells is actually much more complex than that 
simply described by diffusion, it is affected by chemotaxis and haptotaxis. The dynamics 
of tumor-host interactions is determined by many complex cellular and extracellular pro- 
cesses, which include normal epithelial and mesenchymal cells as well as the extracellular 
matrix in addition to the immune response. This interaction is very complex and highly 
variable. Normal cells adjacent to the tumor are often induced to promote tumor growth 
by releasing proteolytic enzymes to break down the extracellular matrix or release growth 
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factors that enhance tumor prohferation. In addition, it appears that tumor cells may under 
some circumstances transform into mesenchymal cells producing new populations of rela- 
tively normal appearing cells that support tumor growth. Tumors are extensively infiltrated 
by immune cells which may constitute as much as one third of its volume. Both the tumor 
phenotype and the tumor environment are very heterogeneous. The former is the result of 
accumulating random mutations, variable enviromental selection forces and perhaps restric- 
tion of proliferate capacity in non-stem cell components of the tumor. In addition, the tumor 
environment is extremely heterogeneous primarily due to disordered angiogenesis and blood 
flow. These facts underline the huge complexity of the problem, and remind us that the 
equations appearing below do not constitute a fundamental description of tumor growth, 
but a statistical approach to some of its properties. 

The continuum equation which describes the MBE universality class (also known as the 
Mullins- Herring equation is the following SPDE 

dth = -KV^h + F + r]{^,t), (1) 

where h is the interface height, K is the surface diffusion coefficient, and ri{'x, t) is a gaussian 
noise with zero mean and correlations given by < ri{'x,t)r]{'x' ,t') >= D6{x — x.')6{t — t'). 
The term F has the dimensions of a velocity, and in this case should be interpreted as the 
product of the mean cell radius and the cell division rate. The critical exponents can be 
extracted from this equation simply by power counting. If we neglect for a moment the 
velocity F, and we perform the transformations x — bx, t b^t, and h — * b'^h in the one- 
dimensional case, we see that the only values of z and a that yield scale invariance are z = 4 
and a = 3/2. Since in this case a > 1, the system is super-rough and it is characterized by 
the set of critical exponents: a = 3/2, aioc = 1, z = 4, j3 = 3/8, and (3* = 1/8 0|, which 
were found to be compatible with those measured in experiments on tumor growth js, Ql- 
Another additional nice feature of the Mullins-Herring equation is its simplicity: since it is 
linear, it can be exactly solved by means of a Fourier transformation. Another properties of 
this equation, however, make it not so suitable for describing tumor growth. It describes the 
growth of a surface from a planar substrate of fixed size, while actual tumors show radial 
growth with their size continuously increasing in time. It is thus important to derive a SPDE 
able to describe MBE physics with the correct geometry and spatiotemporal properties of 
tumor growth. 
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In order to get the correct theoretical description of tumoral growth we will need to borrow 
some elements from differential geometry. On the other hand, these geometrical concepts are 
common in the formulation of stochastic growth equations in reparametrization invariance 
form 9]. The equation of growth of a general riemannian surface reads 

dtf{s, t) = n(s, t)r[f(s, t)] + $(s, t), (2) 

where the d+1 dimensional surface vector r(s, t) = {rQ,(s, t)}^jl\ runs over the surface as 
s = varies in a parameter space (in the following, latin indices vary from 1 to ci and 

greek indices from 1 to + 1). In this equation n stands for the unitary vector normal at 
the surface at r, F contains a deterministic growth mechanism that causes growth along the 
normal n to the surface, and $ is a random force acting on the surface. In our case the 
deterministic part should include a term modelling cell diffusion in the tumor border. When 
surface diffusion occurs to minimize the surface area the corresponding term in the equation 
is Q: 

r, = -KAblH, (3) 
where Abl the Beltrami-Laplace operator 

Abl = -^d.{V99''d,), (4) 

Qij is the metric tensor and g is its determinant, di = d/ds'^ is a covariant derivative, and 
H = n ■ A^^r is the mean curvature. Summation over repeated indices is always assumed 
along this work. Finally, the unitary normal vector is given hy h = g^^^'^difx ■ ■ ■ x S^r. In 
the case of the (l+l)-dimensional Monge form (or what is the same, the parametrization 
corresponding to a planar substrate) we have r = (x, h{x)), the unitary normal vector takes 
this times the form 

the metric tensor is given by (note that for this particular case the metric tensor is a scalar) 
^ = 1 + {dxhY. Thus, the resulting mean curvature is 

^ dlh 

[i + (a^/,)2]3/2- y) 

The corresponding contribution to the drift reads 

^ 3[-i + b{dxhf]{dlhf - iQ[d^h + {dxhf]dlhdlh + [1 + {d^hffdth 

[1 + (9,/i)2]9/2 ' y'> 
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that constitutes an expression far more complex than that of the Mulhns-Herring equation. 
However, we can hnearize this expression about the derivatives of h to get Tg = —Kd^h, 
recovering the famihar drift of the MBE equation. This is the so called small gradient 
expansion, and is valid if sharp changes in the interface are absent [9]. The other contribution 
to the dynamics comes from particle input Tp = n- < F >, where F is the flux of cells 
generated at the interface. If we suppose that cell generation is an isotropic process we find 

< F >= Fh, which implies Tp = F. However, the input of new cells F is a random process, 
giving rise to a stochastic contribution to the dynamics under the form of a noise 77 = n ■ $, 
where $ = F— < F >. In summary, this implies that the stochastic force fulfills ^ = nr) 
and < rj >= 0. Rearranging all the terms we recover the one-dimensional Mulhns-Herring 
equation for MBE growth: 

+ F + (8) 

where the noise, r]{x,t), has been assumed to be gaussian with correlation given by 

< r]{x,t)ri{x' ,t') >= D6{x — x')6(t — t'). The drift of this equation comes originally from 
Eq.Q, which expresses the "diffusion of the mean curvature" of the surface. This corre- 
sponds, also, to a homogeneization of the pressure; let us show this fact as follows. Normal 
cells adjacent to the tumor exert force against new born tumoral cells. Concavities are 
sourrended by a higher number of normal cells than convexities, and thus feel more pres- 
sure. Tumoral cells move along the tumor edge driven by the surface forces, what causes 
the effect of a diffusion: tumoral cells are redistributed from convexities to concavities (this 
fact will be shown explicitly below by means of linear stability analysis). The "equilibrium" 
distribution corresponds to the spherically symmetric form, which presumably implies the 
homogeneization of the pressure all along the tumor edge; this same form implies the min- 
imization of the mean curvature of the surface. The specific form of the cuartic derivative 

in this equation Ms been deduced p>_o.og,caUyBa 

Now that we have identified the physical mechanisms that have led us to Eq.(jH)), we are 

in position to derive SPDEs describing the same physics but with geometrical properties 
compatible with those of a tumor. For the case of (l+l)-dimensional circular model in polar 
coordinates we have r = {r{9) cos{6), r{6) sin(^)), the unitary normal vector reads 



the corresponding metric tensor is given this time by ^ = r^ + (96ir)^, and the mean curvature 
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We are now ready to derive the diffusive drift 



+ 7^ h (11) 
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correspondingly hnearized with respect to the different derivatives of r{6). The term con- 
taining the second derivative of r is irrelevant in the renormalization group sense, so we can 
neglect it to obtain 

K d'^r 

r. - (12) 

The stochastic term comes from the force in the same way as in the last case $ = hrj. 
The noise 77 is a gaussian variable with zero mean and correlation given by 

< 7]{s, t)r]{s', t') >= n„(s, t)np{s, t)D''^^—^6{t - t'). (13) 

V 9 

We now can write the SPDE for tumor growth in (l+l)-dimensions 

where the noise ri{6,t) is gaussian, with zero mean, and correlation given by < 
T]{6, t)r]{6', t') >= D5{9 — 9')5{t — t'). As indicated above, we have assumed that cell genera- 
tion is isotropic, which implies = DS""^ . It is important to note that this time the noise 
is multiplicative, and that it must be interpreted according to Ito, since all the deterministic 
contributions to the drift have been already extracted 9]. Another desirable characteristic 
of this equation is that the variable 9 only varies in [0, 27r] at any time, which represents an 
advantage with respect to using a different coordinate, as for instance the arc length. The 
arc length (a magnitude more similar to x in Eq.(jH))) vary in an interval which depends on 
time (because the tumor grows), what makes more difficult to study the scaling properties 
of the model. We can determine the critical exponents by power counting. The arc lenght 
of a circumference is / = r9, and taking into account that it scales as / — > hi, we deduce that 
the angle scales as ^ ^ b^^°'9. The other two variables scale as r — > 6"r and t —>■ b^t; direct 
substitution reveals that Eq. p4|l is in the MBE universality class. Of course, the range of 
validity of this equation assumes that the interface shows neliglible overhangs in the radial 
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direction compared with the size of the tumor, a fact that has been observed experimentally 
in many cases P]. 

Now we deal with a nonlinear equation, in contrast to Eq.®, that cannot be solved 
simply by means of a Fourier transformation. We can instead employ different techniques 
in order to get some insight into it. For a big enough tumor we can approximate the first 
moment of the radius by the solution the mean-field version of Eq.()14j). i.e., neglecting the 
noise term. On the other hand, the condition of possibility for formulating a continuous 
equation as an adequate description of a tumor is that it is composed by a sufficiently high 
number of cells. So describing the tumor with a continuous equation is in the same order 
of approximation of considering the mean-field level for computing the first moment of the 
radius. The first step in the analysis is to note that the deterministic version of Eq. (jl4p 
admits radially symmetric solutions of the form r{6,t) = R{t) = Ft + Rq, where Rq is the 
radially symmetric initial condition. It is easy to show the linear stability of this solution by 
substituting r{9,t) = R(t) + p{9,t) in Eq. (fT4|l with D = 0, where p is a small perturbation. 
The resulting equation for p is 

dt {Ft + RoYde^' ^ ' 

Since the function p is 27r-periodic in the Q variable we can express it exactly in terms of a 
Fourier series 

oo 

p(Q,t)= PnWe^"^ (16) 

n=— oo 

and by direct substitution in Eq. (fT3|) we see that the Fourier modes obey 

dpn —Kn* 



dt {Ft + Ro) 
We can integrate exactly this equation to obtain 

-Kn'^ 

Pn{t) = pn{to)exp ' 

where 



rPn. (17) 
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{Ro + FtoY {Ro + Ft) 



(18) 
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Pn{to) = ^l pieMe-'^^'de, (19) 



and thus we see that the perturbation decreases in time provided t > to [l(f\. It is also 
important to note that Eq.()18|) might be interpreted as the response of the tumor to an 
external perturbation. Stochastic generation of new cells drives the tumor away from the 



radially symmetric form, while surface diffusion tries to restore it; redistribution of cell 
density after radial symmetry breaking follows the law Eq. lfTHj) . 

Our next step will be to derive the corresponding equation for the growth of a (2+1)- 
dimensional interface. We can parametrize the two-dimensional surface by means of the 
vector f = {r{6, 0) sin(^) cos(0), r(^, 0) sin(^) sin(0), r(^, 0) cos(^)). Implying that the metric 
tensor is 

Border sin^(6') + (S^r) 
The mean curvature can be derived from the metric tensor, however, the expression is 
so cumbersome that it cannot be handled with simplicity. We can instead linearize this 
expression about the different derivatives of r to get 

dlr 
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(20) 



H = — 



(21) 



Collecting the results we see that for the (2+l)-dimensional model for spherical growth the 
SPDE reads 



dr 
di 



+ 



+ 



+ F + 



1 



r](e,4>,t), (22) 



r\/ I sm( 

where the noise T]{6,(j),t) is gaussian, with zero mean, and correlation given by < 
ri{9, 0, t)rj{6' , 0', t') >= D5{6 — 6')5{(t) — (t)')S(t — t'), and it must be again interpreted accord- 
ing to Ito. As happened with the (l+l)-model, we get the desirable characteristic that the 
variables ^, G [0,27r]. In this case we see again that there exists the radially symmetric 
solution for the deterministic version of Eq. ()22|l : r{6, 0, t) = R{t) = Ft + Ro, where Ro is the 
radially symmetric initial condition. We can analize its linear stability by substituting the 
solution r{6, 0, t) = R{t) + p{6, 0, t) in Eq. ()22|) with D = 0, where p is a small perturbation 
27r-periodic in both variables 6 and 0, and thus can be represented in the form of a Fourier 
series 

oo 

p{9,M= Yl Pn,™(t)e*"'+'"^'^. (23) 



n,m=—oo 



The Fourier modes obey the ordinary differential equation 

dpn,m _ -K 

dt ~ {Ft + RoY 
that can be integrated to yield 



H — m^n^ H — m'^ 
3 3 



Pn,r 



Pn,m(t) = Pn,m{to) exp 



-K 



3F 



1 



(i?o + Ftor {Ro + Fty 



(24) 



(25) 



where 

f27r r'iir 



PnAh) = ^l / p{e,<i),t)e-'^'-''^'t>ded<i), (26) 
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implying the stabihty of the solution provided t > to [ll|. This is, as in the former case, 
the law that imposes density redistribution after radial symmetry breaking due to stochastic 
generation of new cells. 

In conclusion, we have derived the equations of growth of the (1+1)- and (2+1)- 
dimensional tumor interfaces containing the physics of MBE in the appropiate geometry. 
These equations provide us a description of the tumor in a coordinate system that is polar, 
and since all the coordinates are angles we have the additional advantage that they vary in 
intervals that are independent of time. They correctly predict the constant velocity growth 
regime found experimentally during the initial phase of growth, and a linear stability anal- 
ysis of radial solutions allowed us to quantitatively estimate the law of density distribution 
of new generated cells. However, latter stages of growth are characterized by certain dece- 
laration of the growth rate; this fact is not captured by the present model, and it will be 
studied in the future. We have assumed all along this work that the tumor is composed of 
an enough large number of cells so that the hydrodynamic description by means of contin- 
uous equations makes sense. If we want to describe small tumors a kinetic approach to the 
problem becomes necessary, as for instance a master equation formulation. Master equation 
descriptions of growth models are already present in the literature [3] , and may be adapted 
for the present case of a tumor. Furthermore, we can project the master equation for the 
cell population into a SPDE via field theoretic arguments in order to recover a more 
similar theoretical approach to the one presented here. These and other questions will be 
the object of future research. 
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